针对弹簧-阻尼(spring-damp, SD)并联模型描述机器人接触作业时需要耗时调节刚度、阻尼系数, 并在接触末端额外安装力传感器等缺陷, 基于互补问题, 描述空间机械臂末端与目标的单边接触, 推导具有紧凑数学形式的空间机器人接触动力学模型。采用 Lemke 算法设计动力学模型的数值计算方法, 并通过数学仿真验证动力学模型的有效性。